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ABSTRACT 

We present a simple, redshift-independent analytic model that explains the local Fun- 
damental Metallicity Relation (FMR), taking into account the physical processes of 
star formation, inflow of metal-poor intergalactic medium (IGM) gas, and the outflow 
of metal rich interstellar medium (ISM) gas. We show that the physics of the FMR 
can be summarised as follows: for massive galaxies with stellar mass M* ^ 10^ ^M©, 
ISM metal enrichment due to star formation is compensated by inflow of metal poor 
IGM gas, leading to a constant value of the gas metallicity with star formation rate 
(SFR); outflows are rendered negligible as a result of the large potential wells of these 
galaxies. On the other hand, as a result of their smaller SFR, less massive galaxies 
produce less heavy elements that are also more efficiently ejected due to their shallow 
potential wells; as a result, for a given M*, the gas metallicity decreases with SFR. 
For such galaxies, the outflow efficiency determines both the slope, and the knee of 
the metallicity-SFR relation. The model is then successfully matched to results from 
numerical simulations including metal enrichment and feedback at higher redshifts. 

Key words: galaxies:high-redshift - fundamental parameters - evolution - abundances 
- stellar content 



1 INTRODUCTION 

Observationally, the stellar mass and the gas phase metallic- 
ity are two of the most fundamental properties of galaxies, 
reflecting the amount of baryons locked up in stars, and the 
amount of metals present in the interstellar medium (ISM) 
gas, respectively. They also provide insights into a number 
of parameters pertinent for understanding galaxy evolution 
including the stellar initial mass function (IMF) which de- 
termines the amount of metals produced and returned by 
stars to the ISM, the star formation history which deter- 
mines the total stellar mass and metal content, and the 
amount of metal rich gas ejected/mass of metal poor gas ac- 
creted. The existence of a mass-metallicity relation has now 
bee n observationa ll y established by man y authors includ- 
ing iGarnettl (l2002l'l.lTremonti et a"jTll2004l '). [Cowie fc Barged 
l|2008ll and IPerez-Montero et al.l (|200 9l'). The origin of this 
relation, however, has been widely debated and a number of 
different explanations have been proposed fo r its existence , 
including: the ejec tion of metal rich gas (e.g. lGarnettll2002l : 
[s^ oni et al]|20ld ) , a depende nce of the star form ation rate 
(SFR) on the galaxy mass (e.g. Calu ra et al1l2009l ). infall of 
met al poor intergalactic medium (IGM) gas into the galaxy 
(e.g. iDekel et al.l2009l : lMannucci et al.ll2009l ). and a balance 
between ISM metal enrichment due to star formation and di- 
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lution due to infall of metal poor IGM gas (|Finlator fc Davi 
l2008h . 

Recently, u s ing o bservations of galaxies at < 2 < 4, 
iMannucci et al.l (|2010l ) have shown that the mass-metallicity 
relation arises as a consequence of a more fundamental 3D 
relation, referred to as the fundamental metallicity relation 
(FMR), relating the stellar mass, gas-phase metallicity and 
the SFR; indeed, if infall is the main driver for star forma- 
tion, and it is star formation that drives outflows, such a 
relation is only to be expected. 

In this work, our aim is to build a simple and yet phys- 
ically consistent analytic model that explains the observed 
FMR, taking into account the relevant physical processes 
of gas inflow/outflow, star formation and metal production. 
The main advantage of the model presented here lies in the 
fact that its results are essentially independent of the red- 
shift considered for 2^6, for reasons that are explained in 
what follows. We also show that the results from this ana- 
lytic model are in excellent agreement with those obtained 
from state-of-the-art cosmological simulations. 



2 BASIC PHYSICAL INTERPRETATION 

As explained above, our first aim is to build the sim- 
plest model that can account for the observed FMR trends. 
The observable quantities determining the FMR of a given 
galaxy, i.e. stellar mass (M,), SFR (xp), and oxygen abun- 
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dance (often used as a proxy for metallicity) , X — Mo/Mg, 
where Mo and Mg are the galaxy oxygen and gas mass, 
respectively, are related by a simple set of evolutionary or- 
dinary differential equations: 
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assumed to be proportional to the gas mass. The two con- 
stants, R and y, represent the returned fraction from stars 
and the the yield per stellar generation, respectively, and 
are dependent on the IMF; whenever numerical values are 
required, we use {R,y) = (0.79,0.0871) consistent with a 
Salpeter IMF for a lower (upper) mass limit of 1 (100) 
Mq. We assume that both the outflow, 'w{A'I)tp, and in- 
fall, a{M)il>, rates are proportional to the SFR, where M 
is the total (dark-l-baryonic) galaxy mass; while the former 
assumption is easy to understand, the latter has been made 
to obtain the simplest solution to the equations described 
above. We further pose that the metallicity of the infalling 
gas is X. Simple algebraic manipulation of the above 

equations gives the general solution 
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where — Mg/Mgo = 
be expressed as Mgo 

and a = a[R — 1 + a — We now interpret the physical 

meaning of the above relation. 

Firstly, we note that reproducing the observed decreas- 
ing trend of X(^) for the smallest galaxies requires a < Q, 
i.e. the condition a < w + {1 ~ R) has to be satisfied. 
Note that this condition also ensures that ^ 1 at any 
time. In addition, if /i <C 1, as is the case for large galax- 
ies in which most of the gas has been turned into stars, 
X{'ip) ~ y{l — R)a^^ = constant, as indeed observed. This 
property can then be used to normalize the value of a for 
the largest galaxies to the observed oxygen abundance, X°''" 
by requiring a = y{l — R)/X°'"'. Physically, the larger the 
value of a, smaller is the value of X, i.e. the metallicity of 
the galaxy is diluted by low metallicity infalling gas. We note 
that the no-infall case (a — 0) would not admit solutions in 
which metallicity is independent of the SFR; such constancy 
implies that some infall is required. 

Outfiows affect the solution through the term /i"" in 
two ways: (a) by setting the position, -;/)*, of the knee in the 
X('i/)) curve, above which the metallicity starts to drop at 
a fixed M, ; this is determined by the condition fi = 0.5, or 
ip* = e*(l -|- U! — a)A/*. Thus, for large outfiow efiiciencies 
(w), the knee shifts to lower values of ^; at fixed w, less 
massive galaxies deviate from a constant X{tl)) relation at 
lower i/) values, (b) The outflow efficiency w also regulates 
the slope of the curve through the power a: as w is increased 
the slope becomes steeper. The combined effect is that, for 
a given value of ^p, low-mass galaxies are more metal poor 
than massive ones. 

By fltting t he X j^p) curves for various A/* to the 
iMannucci et al.l (|201Cll ) data, we can also tune the depen- 
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Figure 1. Oxygen abundance versus SFR for galaxies of dif- 
ferent stellar masses, in tlie range M, = IO^ '^^+'^'^^-'Mq, with 
j = 0, .., 14, from the bottom to the top curves; the shaded re- 
gion shows the ox ygen abundance-S FR relation inferred obser- 
vationally by Ma nnucci et al ] l|2010l 'l for local SDSS galaxies for 
the same stellar mass range, i.e 1O9-25-"-35M0. Points show the 
evolutionary tracks of two representative galaxies with initial gas 
masses of IO^'^'^^^Mq respectively, and numbers along the track 
indicate the values of r = e»t at selected ages. 



deuce of the only two free parameters of the model: a(M) 
and w{M). In a previous paper (see Fig. 6, iDaval et al.l 
,2009), using the same numerical simulation, we have shown 
that Mt (X M; as a result, we can assume that the two quan- 
tities track each other. We find that the best fit to the data 
require 

Ina = -0.43 



O.O51n(Af./lO^'''^")M0, 



Inw = 1.76 - O.331n(Af./lO^ °)M0, 



(5) 
(6) 



Hence, the accretion efficiency is essentially independent of 
the galaxy mass, implying that the accretion rate scales al- 
most linearly with M. The ejection efficiency is small in large 
galaxies {w « 1.0 for Af« — IO^^'^^Mq) and rapidly increases 
with decreasing Af, , reaching w = 4.8 for the smallest ones 
(A/, — 10^'^^Mq). It is interesting to note that the mass 
dependence we find, w oc M^^^^ cx v^^ , is the one expected 
if the outflow is momentum-driven. 

The metallicity-SFR relation (solid curves) as a function 
of the galaxy stellar mass is shown in Fig. [1] along with the 
evolutionary tracks (points) of two representative galaxies 
whose initial gas mass is Mgo = 1O^°'^^M0 (time runs from 
right to left along the open circles). Galaxies start from low 
O/H, low stellar mass, and high SFR values and then move 
towards higher metallicities as their stellar mass increases. 
However, depending on their mass they move at different 
velocities along the track: the most massive observed galax- 
ies are very evolved with a speciflc age (i.e. in units of the 
star formation timescale eZ^ = 0.62 Gyr) r — eU ^ 4, and 
move essentially along constant metallicity curves; smaller 
objects are younger (1 < r < 2) and are still gently building 
up their metal content. 
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2.1 Remarks 

The simple model presented above reproduces the locally 
observed FMR extremely well, essentially through the sin- 
gle, redshift-independent Eq. 3] Time (or z) only specifies 
when in the course of its evolution a galaxy reaches a given 
(X, ijj) point in the plane. As we see later, the same X ~ tp 
trend emerges from numerical simulations at high-z, thus 
supporting the idea that the FRM is indeed a fundamen- 
tal property of galaxies due to the local interplay of star 
formation, infall and outflows. 

Although t hese c onclusions are similar to those reached 
bv lDave et all l|201ll . Dll), there are important conceptual 
differences. Firstly, Dll make the Ansatz that the galaxy 
system is always in a stationary state, i.e. the infall rate ex- 
actly balances the gas lost in outflows and star formation; 
this is equivalent to setting dMg/dt = in Eq. [2] above, or 
a = (1 — -R) -I- w. Analogously, dX/dt = is imposed on the 
metallicity (Eq. [3}, yielding X ^ y{l ~ R)a~^ in Dll. By 
inspecting Eq. |4l we see that this expression approximates 
the relation for large galaxies with <C 1, but fails to catch 
the metallicity decrease with increasing SFR seen in smaller, 
outflow-dominated galaxies. A posteriori, the values for the 
infall and outflow rate we find by fitting the 2 = FMR 
data are such that a < (1 — R) + w, due to the large out- 
flow rates deduced. Moreover, Dll impose a given external 
infall rate derived from numerical simulations (devised to 
model relatively large galaxies), while we do not make any 
assumption on its mass and redshift dependence; in both 
models though, the infall rate is proportional to t/;. It is not 
clear at this stage how the Dll model would perform when 
compared to the FRM data reproduced here, but it would 
be useful to compare the two approaches in detail. As a final 
point, we stress that our set of time-dependent ODEs (2)- 
(3) require an initial condition Mgo 7^ in order to lead to 
physically meaningful solutions. This implies that our SFR 
always exponentially decrease in time, starting from a finite 
value. There could instead be a relatively short (< 1 Gyr) 
transient period corresponding to the build-up phase with 
increasing SFR (Dll), which is not captured by our simple 
assumptions. Thus a more complex version of the model (or 
numerical simulations, see below) should be used if one is 
interested in the FMR at z ^ 6. 



3 SIMULATING THE FMR 

As argued above, the analytic model presented in Sec. [5] is 
expected to work well for galaxies that have already assem- 
bled their gas mass, once the transient gas build-up phase 
lasting about a Gyr is complete. We now check the valid- 
ity of such a statement using a state-of-the-art cosmological 
simulation with a periodic box size of 75h~^ comoving Mpc. 
This is n ow briefly described an d interested readers are re- 
ferred to lTornatore et all (|2010l ) for a complete description: 
the simulation has been carried o ut using the TreePM-SPIf 
code GADGET-2 ('Springell bOOSi ') with the imp lementation 
of chemodynamics as described in iTornatore e t al. {200"^. 
The adopted cosmological model corresponds to a ACDM 
Unive rse, consistent with th e 5-year analysis of the WMAP 
data (|Komatsu et al.ll2009l ). The run includes metallicity - 
dependent radiative cooling (jSutherland fc Dopital Il993h , 



has an effective model to describe star formation from a 
multi-phase ISM and a prescription f or galactic winds trig- 
gered by supernova (SN) explosions (jSpringel fc Hernguistl 
,2(203). The IMF is taken to be Salpeter in the range 1 — 
100 M0 . Metals and energy are released by stars of different 
masses by pro perly accounting for mass-dep endent lifetimes 
as proposed bv lPadovani fc Matteuccil (119931'). The code uses 
the metallicity-dependent yields from Wooslev fc Weaveil 
(19951 : the yields for SNIa and asymptotic giant branch 
stars have been taken from Ivan den Hoek fc GroenewegenI 
( jl997l ). For each galaxy, i dentified as a grav itationally bound 
group of particles (see ISaro et all l2006l l. in each of the 
snapshots of interest, we obtain the total halo/gas/stellar 
mass {Mh/Mg/Mt), the SFR (1/)) and the mass weighted 
gas/stellar metallicity {Zg/Z^); z f» 4.5 is the lowest red- 
shift snapshot available due to simulation constraints. When 
compared to the standard Sheth-Tormen mass function 
(|Sheth fc TormenI Il999l '). the simulated mass function is 
complete for halo masses Mu 1? 10^'' M© in the entire red- 
shift range of interest; galaxies above this mass cut-off are 
then referred to as the "complete sample" and used for all 
the calculations presented in this paper. 

Comparing the FMR inferred from the numerical sim- 
ulation at z « 4.5 to the analytic one presented in Sec. 
(2] we find that these are in excellent agreement for M* « 
■j^qS. 5-10.45 j^^^ as shown in Panel (a) of Fig. [1 the analytic 
resuhs for three specific mass bins, M* = io«-45.9-25,io.3^^ 
have been plotted in the same panel to guide the eye. 
Such agreement lends support to a scenario wherein galaxies 
evolve from a short-lasting transient gas-accretion phase to 
one in which an interplay between SFR, inflows and outflows 
are the dominant physical processes shaping the FMR: as ex- 
plained in Sec. [2] outflows affect the metallicity-SFR relation 
such that for a given value of M, , galaxies with the largest 
SFR are the most metal poor and the metallicity-SFR re- 
lation flattens with increasing A/* due to a decrease in the 
outflow efficiency; however, the fiattening of the FMR is not 
possible without at least some inflow of metal-poor gas. 

The same FMR trends hold also at z f» 5.7, 6.6 as 
seen from panels (b,c) of Fig. [21 albeit that as expected 
in a hierarchical structure formation scenario where large 
galaxies build up from the merger of smaller systems, the 
stellar mass range narrows with increasing redshift, from 
M, lO*-^-^'' *'M0 at z « 4.5 to « lO^'^-^ '^Mo 
at z ~ 6.6; as a consequence of the narrower M, (and 
SFR) range, the gas metallicity range also decreases to 
^ 0.05 - 0.6(0.06 - 0.4)Zq at z « 4.5(6.6). Finally, al- 
though our simulations reproduce the observed metallicity- 
SFR trend extremely well at all the redshifts considered, we 
do not see a fiattening of such relation, as seen in the z = 
sample (Fig. 1). This is readily explicable considering that 
the fiattening occurs for M* ^ lO^'^M©, independent of red- 
shift (see also Sec. [2]); such galaxies have not yet assembled 
by the lowest redshift presented here, z « 4.5, in the volume 
sampled by the simulation. 

As expected from the discussion in Sec. (2] for a given 
SFR, galaxies with larger M« have a larger gas metallic- 
ity by virtue of their larger potential wells that prevent 
metal loss due to outfiows; this trend holds for all redshifts, 
z ~ 4.5 — 6.6. Again, the SFR range reduces with increas- 
ing redshift, as a consequence of the narrowing of the stellar 
mass range such that ip ~ 0.6 — 200MQyr~^ at z « 4.5 and 
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Figure 2. As a function of the SFR (i/i), we show the gas metallicity, 12 + log(0/H), in different stellar mass bins (M,), as obtained 
from numerical simulations. Panels from left to right show the simulated results for z 4.5 — 6.6, as marked. Filled squares in the 
left-most panel show the SFR-gas metallicity relation obtained using the analytic model presented in Sec. [2] for three stellar mass bins, 
Log(M, ) = 8.5,9.25, 10.3 Mq from bottom to top respectively. This shows the simulated and analytic results are in excellent agreement 
at z ~ 4.5. See text for details. 
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Figure 3. As a function of the stellar mass (M*), wc show the gas metallicity in units of 12 + log(0/H), in different SFR {if}) bins, as 
obtained from numerical simulations. Panels from left to right show the simulated results for z Ri 4.5 — 6.6, as marked. 




reduces to ^ ^ 1-6 — 4OM0yr~^ at z f» 6.6. Finally, due 
to a lack of galaxies with M* 10^^ M0, again, we do not 
see a clear flattening of the metallicity-stell ar mass relation, 
as obs erved for such galaxies at z = by iMannucci et al.l 
(|20ld ). 

Finally, we show the gas metallicity as a function of the 
param eter = log(M,) — 0.321og(i/)) that IMannucci et al.l 
l|2010l l have introduced to minimize the metallicity scat- 
ter. We find that the simulated galaxies follow a tight 
metallicity-(^ relation for all the redshifts considered, as 
shown in Fig. [4] When compared to the relation inferred 
using local SDSS galaxies, z « 4.5 — 6.6 galaxies essen- 
tially follow the same slope on the metallicity-^ plane, al- 
though as expected, both the metallicity and (" ranges de- 
crease with increasing redshift; the metallicity scatter be- 
tween z ~ 4.5 — 6.6 is about 0.35 dex at most (at ( « 9.0). 
Further, at 9.7, the simulated metallicit y at z ~ 4.5 falls 
short by about 0.3 dex of that measured bv lMannucci et al] 
for local S DSS galaxies. This is in slight tension with 
IMannucci et al.l (|2010l ). who find a metallicity difference of 



about 0.6 dex at ~ 9.7 between z ~ 3 and local galaxies, 
although they find no evolution in the metallicity with ^, 
for z — — 2.5. Possible reasons for such a discrepancy are 
discussed briefly below in Sec. |31 



4 CONCLUSIONS 

We have presented a very simple model that provides 
a straightforward explanation of the 3D FMR relating 
the stell ar mass, SF R and gas metallicity observed by 
Mannucc i et al.l l|2010t ) in the local Universe. Starting from 
a galaxy that has already assembled its gas mass after a 
transient gas-accretion phase, assuming the metallicity of 
the infalling IGM gas to be much lower than that of the 
ISM gas, and assuming both the inflow and outflow rates 
to be proportional to the SFR, we calibrate the only two 
model free parameters (the inflow and outflow rates) to the 
local FMR data. Using this model, we show that both in- 
flows and outflows are required to explain the observations, 
although their importance is mass dependent: since we as- 
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Figure 4. As a function of tlie parameter <^ = log(M«) — 
0.32 log('!/'), we show tiie gas metallicity, 12 + Iog(0/H), for 
z ~ 4.5 — 6.6, as obtained from numerical simulations ; the solid 
line shows the measurements of iMannucci et ah I ll2010l'l for local 
SDSS galaxies. 



made. Finally, at = 9.7, we find a smaller metallicity differ- 
ence (about 0.3 dex) between the z ~ 4.5 simulated results 
and z = Q observa tions, c ompa red to the value of 0.6 dex 
found by.Ma nrm c ci et al.l (|2010l ) between the data a,t z — Q 
and 3. Although such a difference could arise from observa- 
tional biases at z « 3, i.e. small number statistics (16 ob- 
served galaxies), or the use of the OIII line t hat is biased to- 
ward s lower metallicity galaxies (see Sec. S.l lMannucci et al] 
|2010D . it could also be in part attributed to the difference 
in the simulated and observed mass ranges, and the effects 
of the resolution and feedback implementation in numeri- 
cal simulations, on the SFR and gas metallicity. However, a 
much larger observed data set is required before any such 
meaningful comparison can be made. 
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sume the inflow rate to be proportional to the SFR (and 
hence M,), galaxies with the largest M* accrete the largest 
amount of metal-poor IGM gas, which is balanced by the 
larger amount of metals produced due to star formation; 
this, coupled with the negligible outflow rates due to their 
large potential wells, makes the metallicity-SFR relation es- 
sentially constant for galaxies with A/* ^ lO'^^M©. On the 
other hand, due to their smaller SFR, less massive galaxies 
produce less metals. However, as a result of their smaller 
potential wells, outflows lead to a loss of metal-rich ISM gas 
so that at a given M*, galaxies with the largest SFR are 
the most metal poor; alternatively, for a given SFR, smaller 
galaxies are more metal poor than massive ones. 

Since the analytic model developed is valid for all galax- 
ies that have 'assembled' their gas mass after a transient 
phase of gas-accretion, we can compare its results with 
those from state-of-the-art cosmological simulations; due to 
computational constraints, such comparison is only possi- 
ble down to 2 f» 4.5. However, even at this high-redshift, 
the analytic results are in excellent agreement with those 
from the simulation, providing a validation for our analytic 
model. The same FMR trends hold up to z « 6.6, although 
as expected, galaxies become smaller and have lower gas 
metallicities with increasing redshift. 

We conclude by discussing the main caveats: firstly, as 
already mentioned, we ignore the complication of including 
the initial phase of gas accretion that is important for rela- 
tively large galaxies. However, as is shown by the excellent 
agreement between the analytic model and results from nu- 
merical simulations, such considerations are important at 
high-redshifts (z ^ 5) when galaxies are still building up 
their final mass. Secondly, to find an analytic solution to 
Eqs. 2-4, we have used an infall rate that is proportional 
to the SFR, although we do not make any assumption on 
its galaxy mass and redshift dependence. In the future, we 
aim at exploring models wherein such an assumption is not 
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